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The  objective,  of  this  work  was  to  develop  a  broad  array  of  methods  for  electronic  structure  theory,  with  special  emphasis  or 
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for  molecules  with  unpaired  electrons  has  been  completed  and  is  in  the  most  recent  GAMESS  release.  The  derivation  for 
the  preferred  restricted  open  shell  ZAPT  method  is  complete  and  a  paper  describing  this  method  has  been  published.  The 
corresponding  code  is  now  in  GAMESS.  The  code  for  implementing  Hartree-Fock  with  DDI  has  been  completed,  and  a 
paper  describing  this  code  has  been  published.  We  are  now  turning  to  analogous  developments  for  MCSCF  and  Cl  codes. 

A  scalable  code  for  the  multi-reference  second  order  perturbation  theory  has  been  completed,  a  paper  published,  and  the 
code  is  in  GAMESS.  Other  important  developments  include:  the  derivation  of  gradients  for  multi-reference  second  order 
perturbation  theory,  further  developments  of  our  effective  fragment  potential  (EFP)  method  for  studying  solvation  and  liqui 
behavior,  the  development  of  molecular  dynamics  and  Monte  Carlo  methods  to  facilitate  the  study  of  solvation  and  liquid 
behavior,  the  development  and  implementation  of  a  new  method  for  producing  global  energy  surfaces  from  (cont.  below) 
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3.  Status  of  effort:  A  brief  statement  of  progress  towards  achieving 
the  research  objectives.  (Limit  to  200  words.) 

The  objective  of  this  work  was  to  develop  a  broad  array  of  methods  for  electronic 
structure  theory,  with  special  emphasis  on  improved  parallel  codes.  A  new  highly 
scalable  code  for  frozen  core  second  order  perturbation  theory  (MP2)  gradients  for  closed 
and  open  shell  molecules  has  been  developed  and  is  on-line  for  general  use  at  several 
IBM  SP2  and  Cray  T3E  systems.  The  key  element  of  this  approach  is  our  distributed 
data  interface  (DDI)  that  facilitates  distribution  of  very  large  arrays  across  all  available 
nodes.  This  code  has  been  ported  to  the  cluster  environment,  so  it  can  now  be  used  on 
clusters  of  PCs  or  Macs  running  Linux,  as  well  as  on  clusters  of  high  performance 
workstations.  The  development  of  analogous  codes  for  molecules  with  unpaired  electrons 
has  been  completed  and  is  in  the  most  recent  GAMESS  release.  The  derivation  for  the 
preferred  restricted  open  shell  ZAPT  method  is  complete  and  a  paper  describing  this 
method  has  been  published.  The  corresponding  code  is  now  in  GAMESS.  The  code  for 
implementing  Hartree-Fock  with  DDI  has  been  completed,  and  a  paper  describing  this 
code  has  been  published.  We  are  now  turning  to  analogous  developments  for  MCSCF 
and  Cl  codes.  A  scalable  code  for  the  multi-reference  second  order  perturbation  theory 
has  been  completed,  a  paper  published,  and  the  code  is  in  GAMESS.  Other  important 
developments  include:  the  derivation  of  gradients  for  multi-reference  second  order 
perturbation  theory,  further  developments  of  our  effective  fragment  potential  (EFP) 
method  for  studying  solvation  and  liquid  behavior,  the  development  of  molecular 
dynamics  and  Monte  Carlo  methods  to  facilitate  the  study  of  solvation  and  liquid 
behavior,  the  development  and  implementation  of  a  new  method  for  producing  global 
potential  energy  surfaces  from  sets  of  ab  initio  points,  the  development  and 
implementation  of  both  grid-based  and  gridless  approaches  to  density  functional  theory, 
and  the  development  and  implementation  of  several  MCSCF-based  approaches  to  spin- 
orbit  coupling. 

4.  Accomplishments/New  Findings:  Describe  research  highlights,  their 
significance  to  the  field,  their  relationship  to  the  original  goals,  their  relevance  to 
the  AF's  mission,  and  their  potential  applications  to  AF  and  civilian  technology 
challenges. 

The  objective  of  this  work  is  to  develop  a  broad  array  of  methods  for  electronic  structure 
theory,  with  special  emphasis  on  improved  parallel  codes.  Progress  has  been  made  in 
many  areas. 

Parallel  Developments.  A  highly  scalable  code  for  frozen  core  second  order  perturbation 
theory  gradients  for  closed  shell  molecules  has  been  developed  and  is  on-line  for  general 
use  at  several  MSRC  sites.  This  new  code  uses  the  Distributed  Data  Interface  (DDI)  so 
that  the  large  arrays  do  not  have  to  be  replicated.  This  means  that  much  larger  molecules 


can  be  studied  using  geometries  based  on  correlated  wavefunctions.  The  code  scales  very 
well  for  as  many  as  5 12  nodes  and  has  already  been  applied  to  several  challenging 
compounds.  The  development  of  analogous  codes  for  molecules  with  unpaired  electrons, 
using  unrestricted  wavefunctions,  has  been  completed.  The  corresponding  method  for 
spin  restricted  wavefunctions  has  been  completed.  The  unrestricted  second  order 
perturbation  theory  parallel  code  under  DDI  has  been  completed  and  tested,  and  a  paper 
on  this  has  been  published.  A  paper  on  the  restricted  open-shell  derivation  has  been 
published,  and  the  development  of  the  parallel  code  has  been  completed.  A  general 
overview  paper  on  second  order  perturbation  theory  derivatives  has  been  published.  A 
paper  on  a  DDI  implementation  of  the  SCF  method  has  also  been  published.  Improved 
parallel  MCSCF  codes  have  been  completed,  and  a  distributed  parallel  full  Cl  code  has 
been  implemented  and  a  paper  published.  The  multi-reference  perturbation  theory  code 
(MCQDPT)  has  now  been  implemented  under  DDI.  The  effective  fragment  potential 
(EFP)  code  has  been  made  scalable  under  DDL  In  an  important  related  development, 
parallel  GAMESS  now  runs  on  clusters  of  PCs  and  Macs  (running  Linux)  and  high  end 
workstations.  This  is  more  challenging  than  self-contained  massively  parallel  computers, 
since  the  overhead  due  to  inter-node  communications  is  more  complex.  This  issue  has 
been  solved  to  some  extent  by  using  a  Gigabit  Ethernet  switch  with  large  data  packets. 

We  are  also  exploring  alternative  communications  solutions,  such  as  Myrinet,  SCI,  and 
Infiniband.  These  developments  and  their  applications  have  been  enhanced  by  the 
construction  of  a  new  32-node,  128-CPU  IBM  Power3+  cluster  using  funds  provided  by 
a  DURIP  grant.  As  for  all  other  GAMESS  developments,  we  will  make  our  experiences 
in  developing  scalable  clusters  available  to  all  users. 

Condensed  Phase  Methods.  We  have  already  shown  that  our  effective  fragment  potential 
(EFP)  method  for  solvation  is  excellent  for  water,  in  a  variety  of  applications.  We  are 
now  working  on  extending  the  capabilities  of  the  method  in  several  ways.  We  have 
developed  a  new  approach  for  incorporating  dispersion,  charge  transfer  and  other  higher 
order  terms  into  the  method.  Such  terms  are  particularly  important  for  nonpolar  solvents. 
The  model  has  been  developed  in  such  a  way  that  it  is  equally  applicable  and  accurate  for 
any  solvent.  The  derivation  of  an  expression  for  the  analytic  gradient  for  the  EFP-a6 
initio  interaction  term  is  in  progress.  A  very  important  new  development  (by  co-worker 
Jan  Jensen)  is  a  new  method  for  using  EFPs  across  covalent  bonds.  This  will  facilitate  the 
representation  of  large  substituents,  as  well  as  the  treatment  of  large  biological  molecules. 
We  have  also  interfaced  the  EFP  method  with  two  continuum  methods,  the  simple 
Onsager  reaction  field  and  the  more  sophisticated  polarizable  continuum  model  (PCM). 

As  the  number  of  solvent  molecules  in  the  system  increases,  the  number  of  configurations 
to  be  considered  increases  rapidly,  and  traditional  small  molecule  geometry  optimization 
methods  are  not  effective.  We  are  therefore  developing  both  molecular  dynamics  and 
Monte  Carlo  simulation  codes  so  that  the  configurational  space  can  be  probed  more 
effectively,  not  only  for  minima,  but  for  transition  states  and  reaction  paths  as  well.  This 
is  also  an  important  first  step  in  the  development  of  methods  based  on  the  EFP  model  for 
predicting  bulk  properties  and  super-critical  behavior.  We  are  therefore  incorporating 
Ewald  summations  into  our  MD  method.  Related  to  these  dynamical  methods  is  our 
development  and  implementation  of  a  method  for  converting  the  large  numbers  of  points 


generated  in  ab  initio  trajectories  into  a  global  potential  energy  surface,  using  a  modified 
Shepard  interpolation  approach. 

Major  James  Shoemaker’s  Ph.D  dissertation  in  Engineering  Physics  at  the  Air  Force 
Institute  of  Technology  focused  on  the  development  and  applications  of  an  embedded 
cluster  model  for  treating  surface  chemistry.  A  paper  on  the  theoretical  method,  called 
SIMOMM,  has  been  published  in  the  Journal  of  Physical  Chemistry.  Several  papers 
applying  this  method  have  now  been  published  or  submitted,  and  a  manuscript  that 
describes  applications  to  silicon  carbide  surfaces  is  in  progress.  Extensions  of  this 
method  to  transition  metal  catalysts  are  planned;  however,  molecular  mechanics  is  not 
likely  to  be  a  viable  approach  for  the  bulk,  since  electrons  in  such  systems  are  too 
delocalized.  To  overcome  this  potential  problem,  we  have  begun  to  explore  fast 
multipole  methods  (FMM)  that  scale  linearly  and  are  also  highly  parallel,  so  that  the  bulk 
part  of  the  system  can  be  treated  by  quantum  mechanics.  One  preliminary  paper  has  been 
published  and  the  first  FMM  paper  has  been  published.  The  parallel  efficiency  of  the 
FMM  method  has  been  greatly  improved,  and  gradients  have  recently  been  developed 
and  implemented. 

Other  developments  There  are  many  problems  that  involve  species  for  which  one  cannot 
write  one  simple  Lewis  structure.  Such  species  must  be  treated  with  multi- 
configurational  (MCSCF)  wavefiinctions.  However,  the  computational  demands  for 
MCSCF  calculations  are  very  high.  So,  we  have  developed  a  set  of  methods  for 
converging  such  wavefiinctions  that  range  from  very  rapid  convergence  with  high  CPU 
and  disk  requirements  to  methods  that  converge  more  slowly  but  have  very  modest 
computational  demands.  The  advantage  of  the  latter  is  that  much  larger  molecules  than 
normal  can  be  treated  with  this  method. 

For  species  that  must  be  treated  with  MCSCF  wavefiinctions,  non-traditional  methods  are 
also  needed  to  recover  the  dynamic  correlation.  One  method  that  has  been  implemented 
into  GAMESS  is  multi-reference  second  order  perturbation  theory.  The  gradient 
formulas  for  this  method  have  been  derived,  and  the  implementation  into  GAMESS  is  in 
progress. 

There  are  many  problems  that  involve  multiple  electronic  states  on  multiple  potential 
energy  surfaces.  Since  such  surfaces  can  cross  and  therefore  induce  radiationless 
transitions,  it  is  important  to  have  the  capability  to  predict  non-adiabatic  coupling  among 
these  surfaces.  One  origin  of  non-adiabatic  interactions  among  surfaces  is  spin-orbit 
coupling.  Spin-orbit  coupling  matrix  elements  must  be  calculated  using  multi-reference 
wavefiinctions,  but  until  recently,  there  have  been  no  generally  available  codes  for  the 
computation  of  these  matrix  elements  using  all  electrons  and  any  arbitrary  MCSCF 
wavefunction  for  any  arbitrary  set  of  spin  states.  We  have  now  developed  and 
implemented  such  a  method.  It  will  be  particularly  valuable  for  the  study  of  high  energy 
density  materials  and  for  the  study  of  any  processes  involving  transition  metals,  since 
both  types  of  species  involve  multiple  electronic  states.  Since  such  all-electron 
calculations  are  very  demanding,  two  alternative  approaches  to  this  problem  have  also 
been  develop  and  implemented  into  GAMESS  for  use  with  larger  molecules.  One  of 


these  is  called  the  Zeff  method,  in  which  the  two  electron  term  is  ignored  and  a  fitted 
empirical  parameter  is  introduced  into  the  one-electron  part  of  the  Hamiltonian.  While 
this  is  a  very  efficient  approach,  its  accuracy  is  limited.  However,  there  is  an 
intermediate  method,  in  which  only  those  matrix  elements  connecting  orbitals  in  the 
active  space  with  each  other  are  neglected.  We  have  shown  that  in  this  case  excellent 
agreement  with  the  all  electron  method  is  obtained,  without  any  need  for  fitted 
parameters.  We  have  also  implemented  the  RESC  (relativistic  elimination  of  the  small 
component)  method  for  calculating  spin-orbit  coupling  interactions,  complete  with 
analytic  gradients. 

Density  functional  theory  (DFT)  has  become  extraordinarily  popular  in  recent  years, 
because  it  has  the  same  computational  requirements  as  Hartree-Fock  theory,  but  often  has 
a  level  of  accuracy  that  is  comparable  to  second  order  perturbation  theory.  We  have 
developed  and  implemented  a  gridless  approach  to  DFT,  by  using  the  resolution  of  the 
identity  to  evaluate  the  complicated  integrals.  This  procedure  requires  the  use  of 
auxiliary  basis  sets  to  resolve  the  identity.  After  extensive  analysis  of  the  convergence  of 
several  molecular  properties  as  a  function  of  basis  set,  we  have  recommended  a  series  of 
auxiliary  basis  sets  for  use  with  the  first  row  in  the  periodic  table.  The  next  step  will  be 
to  extend  this  basis  set  development  to  heavier  elements.  We  have  also  implemented  a 
more  traditional  grid-based  DFT  method,  including  all  of  the  most  popular  functionals. 

Another  approach  that  is  very  important  for  large  molecules,  especially  molecules  that 
contain  heavier  elements,  is  the  effective  core  potential  (ECP)  method.  This  method 
replaces  inner  shell  electrons  by  model  potentials,  making  the  calculations  much  more 
efficient.  These  potentials  include  relativistic  corrections  for  the  heavier  elements.  To 
make  these  calculations  more  useful,  we  have  developed  much  more  efficient  methods 
for  obtaining  the  first  and  second  energy  derivatives;  that  is,  the  gradients  and  hessians. 
The  increase  in  speed  is  a  factor  of  5-10,  and  the  method  has  been  implemented  to  take 
advantage  of  parallel  computers. 

As  calculations  on  larger  and  larger  species  are  made  feasible  by  advances  in  parallel 
programs  and  other  innovative  algorithms,  it  becomes  increasingly  important  to  develop 
graphical  methods  for  interpreting  the  results  of  complex  calculations.  This  has  led  to  the 
continuing  development  of  MacMolPlt,  a  graphical  interface  for  GAMESS.  MacMolPlt 
provides  the  ability  to  visualize  (and  rotate  in  real  time)  orbitals,  wavefunctions,  electron 
densities  and  electrostatic  potentials  in  two  and  three  dimensions.  We  are  now  in  the 
process  of  developing  a  Windows  based  version  of  this  code. 

We  have  also  been  exploring  alternatives  and  extensions  to  the  popular  G2  and  G3 
methods  developed  by  Pople  and  co-workers  for  the  accurate  prediction  of  such 
theormodynamic  properties  as  heats  of  formation,  ionization  potentials  and  electron 
affinities.  One  limitation  of  these  methods  is  that  they  are  applicable  only  to  species  that 
are  adequately  described  by  single  configuration  wavefunctions.  This  eliminates  many 
transition  metal  compounds,  diradicals  and  most  transition  states.  We  have  therefore 
initiated,  in  collaboration  with  the  Radom  group,  the  development  of  multi-reference 
analogs  of  the  G2  and  G3  methods.  The  multi-reference  methods  are  based  on  CASSCF 


wavefimctions,  instead  of  Hartree-Fock,  followed  by  multi-reference  perturbation  theory 
(instead  of  MP2)  and  finally  multi-reference  CL  The  initial  paper  describing  this  method 
has  been  published. 
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